Keeping Accurate Time For A Hybrid GPS Receiver And Mobile Phone When Powered Off

ABSTRACT

A hybrid navigation satellite receiver and mobile telephone uses only two crystal oscillators. One that operates a master clock around 27-MHz and that consumes milliwatts of power. The other oscillator consumes only microwatts of power and operates continuously on battery power at about 32-KHz. Only the second, low frequency oscillator is kept running during power “off”. On power “restart”, a real-time-clock counter is consulted to cause an estimate of the GPS system time to be regenerated and supplied to the GPS-DSP to quicken its initialization. The master clock is GPS-calibrated, and the accurate clock is used to drive NCO&#39;s for the mobile telephone part and host CPU.

CROSS-REFERENCE TO RELATED APPLICATIONS

This application is a divisional of U.S. patent application Ser. No. 11/110,200 filed Apr. 20, 2005, which is a continuation-in-part of patent application Ser. No. 10/656,078 filed Sep. 4, 2003, now abandoned. The contents of these applications are incorporated herein by reference in their entirety.

BACKGROUND

1. Field of the Invention

The present invention relates to the problem of slow initializations when time is uncertain, and more particularly to using a GPS-calibrated real-time clock left running during power-down cycles to supply an accurate time estimate for both halves of a hybrid GPS-receiver mobile-phone during initialization.

2. Description of the Prior Art

Most consumer electronic devices are mass produced and their sales are very dependent on how they are priced. One way to offer better prices is to reduce manufacturing costs, e.g., the cost of labor and the components used. Combinations of navigation satellite receivers and communication devices are now available, and many conventional commercial products provide one set of crystals and crystal oscillators for the navigation satellite receiver part, and a separate set for the communication devices.

Global positioning system (GPS) receivers use signals received from several earth-orbiting satellites in a constellation to determine user position and velocity, and other navigational data. A prior art navigation receiver that has just been turned on does not yet know where it is, how much its crystal oscillator is in error, nor what time it is. All these are needed to find and lock onto the satellite transmissions, and so a search must be made of the possibilities. The most probable are searched first to save time.

High-sensitivity GPS receivers are a problem when the initial time or frequency uncertainty is large. Finding signal energy when the signal energy is extremely faint requires making smaller steps and dwelling at each step longer. So having a better initial estimate of the local reference oscillator can improve time-to-first-fix.

GPS receivers with signal levels better that −145 dbm can readily lock onto a strong GPS satellite vehicle (SV) to decode the NAV-data. Such yields the SV ephemeris and position. After than, the total pseudorange needs to be formed from the hardware codephase. Conventional GPS receivers determine the integer millisecond and so-called z-count.

When signal levels are roughly no better than −145 dbm to −150 dbm, a practical high-sensitivity GPS receiver can employ pattern-matching tricks to get a z-count or integer millisecond for an anywhere-fix.

GPS receivers that are locked onto and tracking one or more GPS satellite vehicles (SV's), know time to very high accuracy. This is because the GPS system is based on atomic clocks that set the time and frequency references used. The coarse acquisition (C/A) transmitted by the SV's repeats every one millisecond of propagation wavelength, and so there is a basic ambiguity as to which millisecond a GPS receiver is observing. If the integer millisecond is known, e.g., time is known to better than one millisecond, then the integer ambiguity does not need to be solved. The z-count is known. Skipping the steps to find z-count and set the integer millisecond can save a tremendous amount of time and effort in a GPS receiver working on providing its first navigation solution fix after a cold start.

There is a need for a less expensive-to-produce combination of a navigation satellite receiver and a communication device.

SUMMARY OF THE INVENTION

Briefly, a hybrid navigation satellite receiver and mobile telephone embodiment of the present invention uses only two crystal oscillators. One that operates a master clock around 27-MHz and that consumes milliwatts of power. The other oscillator consumes only microwatts of power and operates continuously on battery power at about 32-KHz. Only the second, low frequency oscillator is kept running during power “off”. On power “restart”, a real-time-clock counter is consulted to cause an estimate of the GPS system time to be regenerated and supplied to the GPS-DSP to quicken its initialization. The master clock is GPS-calibrated, and the accurate clock is used to drive NCO's for the mobile telephone part and host CPU.

An advantage of the present invention is that a system and method are provided for a hybrid navigation receiver and mobile telephone combination to initialize quickly and with only two crystals.

These and other objects and advantages of the present invention will no doubt become obvious to those of ordinary skill in the art after having read the following detailed description of the preferred SPS receivers which are illustrated in the various drawing figures.

IN THE DRAWINGS

FIG. 1 is a functional block diagram of a combination GPS satellite navigation receiver part and communication device system embodiment of the present invention;

FIG. 2 is a functional block diagram of a GPS system embodiment of the present invention that supports the client receiver with information supplied by a network server; and

FIG. 3 is a schematic diagram of a multiple re-use of frequencies in a combination GPS satellite navigation receiver and communication device system embodiment of the present invention.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT

FIG. 1 illustrates a hybrid satellite navigation receiver and communication device embodiment of the present invention, and is referred to herein by the general reference numeral 100. The system 100 has two basic parts, a navigation receiver in the form of a global position system (GPS) receiver 102, and communication device in the form of a mobile telephone 104. The GPS part 102 is either a self-contained positioning model or a measurement model with positioning code executed on a host CPU. Each part can derive accurate system time respectively from the GPS system and the cellular telephone infrastructure. Accurate time is needed to stabilize their respective local oscillators (LOgps and LOcell). Having accurate local oscillators allows each part to find carriers in their respective spectrums without under frequency search times.

Embodiments of the present invention allow a complete shut down of both the GPS receiver 102 and the mobile telephone 104 to save battery power. The GPS receiver 102 comprises a GPS antenna 106 connected to a GPS-RF-IF stage 108. A GPS digital signal processor (GPS-DSP) 110 samples downconverted signals and produces position solutions and highly accurate GPS system time. A 27-MHz crystal 112 controls a master clock 114. The various LO and IF frequencies needed by the GPS receiver are generated from the master clock by GPS oscillator 116. A software controlled crystal oscillator (SCXO) model 118 receives a temperature input from a temperature sensor (T) 120. It also receives accurate GPS time from a timebase 122. The SCXO model 118 is able to compensate the temperature drift of crystal 112. For a complete discussion of this technique, which was incorporated by reference in the parent application and here, see, U.S. Pat. No. 6,509,870 B1, issued Jan. 21, 2003.

During normal operation, the master clock 114 needs milliwatts of power to operate. So this must be shut down when the user inputs an “off” command. Its operation is resumed when the user issues a “restart” command. The problem is, accurate GPS time was not being computed during the off time.

During normal operation, the accurate timebase 122 supplies computed GPS system time to a GPS time regenerator 124. Such also receives an “offset” count from a real time clock (RTC) 126. A low frequency 32-KHz crystal 128 controls the RTC 126, and such operates on microwatts of power. This is three orders of magnitude less power than is needed by master clock 114. The GPS time regenerator 124 is able to compute the error and drift in the low frequency 32-KHz crystal 128. It also receives a count that it aligns with GPS system time. When an “off” command is received, the offset count aligned with the last good GPS system time computed is stored. When a “restart” command is received, the GPS time regenerator 124 adds the current offset count, corrected for error and drift, to the last good GPS system time computed. This results in a regenerated GPS time that is fed back through the accurate timebase 122 to the GPS=−DSP 110. There it is used as a starting point to begin carrier frequency searches and Doppler estimates. It also allows the master clock 112 to immediately begin operation under control of the SCXO model 118. An accurate master clock is therefore immediately available after “restart” to the mobile telephone 104.

Such mobile telephone 104 comprises a communications antenna 130 connected to a communications transceiver 132. A communications DSP 134 decodes and codes mobile voice calls. A host CPU 136 is resident in the portable unit and can be used for other applications. It can also control the exact LO frequencies being synthesized by a first numeric controlled oscillator (NCO1) 138. A second NCO2 140 can be similarly controlled for clocks needed by the CPU 136 or its applications.

FIG. 2 and its associated disclosure was originally incorporated by reference by the parent case, and is repeated here for convenience of the reader. The incorporated material now appears in U.S. Pat. No. 6,714,160, issued Mar. 30, 2004, to one of the present inventors, Paul McBurney. It represents a OMNI-client navigation satellite receiver network, and is referred to herein by the general reference numeral 200. The OMNI-client navigation satellite receiver network includes at least one navigation platform 202 supported by a network server 204.

Each GPS measurement platform 202 typically comprises a GPS antenna 206, a low-noise amplifier (LNA) 208, a GPS surface acoustic wave (SAW) filter 210, a radio frequency (RF) application specific integrated circuit (ASIC) 212 with an intermediate frequency (IF) SAW-filter 214, a digital signal processor (DSP) 216, a reference crystal oscillator 218, and a reference-crystal temperature sensor 220.

An autonomous client 222 can function and provide navigation solutions to a user with next to no help from the server 204. A demi-client 224 needs help, e.g., polynomial models that simplify ephemeris and time bias calculations. A thin-client 226 does not burden its local host with navigation solution processing. It lays off just about all the navigation calculations on the server 204, and basically provides only observational measurements from its point of view of the SV constellation. The navigation solutions are returned for local display if a user is there and wants to see them. In a thin-client 226, the DSP is a shared part with some other non-GPS application. As such, multi-threaded application programs are not needed at the client and only simple program loops are executed.

An OMNI client 227 runs almost fully autonomously, but periodically collects a full set of ephemeredes over the computer network. It further operates during power off to keep its position uncertainty, sigmaPos, under 150-km when it is powered back up. These conditions allow high sensitivity operation where much finer search steps are used to find signal power and each step has a long dwell. The OMNI client 227 also benefits greatly if the crystal oscillator 218 is software compensated using temperature measurements by temperature sensor 220. A real time clock (RTC) is kept running that is accurate to better than one millisecond of true time each time the navigation platform 202 is powered-up.

The local reference oscillator crystal oscillator 218 will have a frequency drift error that will vary as a function of temperature. The reference-crystal temperature sensor 220 is used to measure the temperature of the local reference oscillator crystal oscillator 218. The first use is to collect data to build curve during manufacturing calibration when the navigation platform 202 is initialized and tracking SV's. The subsequent use is to provide an index value so a ninth-order polynomial equation can be computed from stored coefficients while the navigation platform 202 is initializing and trying to lock onto its first SV.

The server 204 typically comprises a number of reference station antennas 228 and 230 that provide GPS signal inputs to a reference station manager 232. A location server 234 can provide support information to the demi-client 224, thin-client 226, and OMNI-client 227 for improved time-to-first-fix and position solution quality. In the case of the OMNI-client 227 operating in high-sensitivity mode, the ephemerid information collected and forwarded by the server 204 enables anywhere fixes with signal levels from SV's under −150 dbm.

A method embodiment of the present invention determines how and when the server 204 is contacted by an OMNI client, e.g., a client and navigation platform 202. The server contact must be infrequent and minimized in many situations because the per-byte communication costs are high or the network is only periodically accessible.

When signal strengths are high, the z-count and BTT are really measured by collecting NAV-data. The BTT is used to clean-up any codephase roll-overs. Generally the sub-20 msec portion should agree. There is a bit more noise on the BTT than z-count. However, the z-count can be off by one millisecond for a short period of time near where the codephase rolls over.

The OMNI client needs a good time source to reduce sigmaTime to under one millisecond. The 50-Hz NAV-data can be used to do a pattern match and indirectly find time. Such can provide a GPS receiver with an adequate time source when a z-count cannot be demodulated. If there is enough confidence in the pattern match, the integer millisecond, intMs, on an SV can also be determined.

If the starting time uncertainty, sigmaTime, is greater than +/−10 msec, the so-called big delta T term (DT) has to be used in the solution fix. Such increases the number of SV's needed by one. A gridFix method can be used when the position uncertainty, sigmaPos, is under 150-km and the intMs on SV's is not available. A no-Z fix type is used when sigmaTime is greater than 10 msec.

A complete GPS almanac, highAccAlm, is sent by the server with ephemeredes rather than almanacs for all GPS SV's. Another complete GPS almanac, mixAccAlm, can be sent by the server and includes older ephemeredes for SV's it is not currently tracking.

Preferably, a WWserver server is implemented that has continuous observability of the complete GPS constellation. It has enough reference stations with adequate spatial separation to view all SV's around the world at the same time. Server 204 represents a local server, LAserver, that has one or more reference stations that observe only a sub-set of the complete GPS SV-constellation. Therefore an LAserver can not provide highAccAlm, only mixAccAlm.

After turning on, the almanac will include ephemeredes that are actually almanacs. After one 12-hour cycle, some of the almanacs will be replaced with ephemeris based almanacs.

The NAV-data from the GPS SV's can be collected directly down to signal levels as low as −145 dbm. Therefore, the ephemeris, z-count and BTT can be derived at this level. SV's at this level can operate independent of the server and also can be used in a fix that has no requirement on starting position accuracy, e.g., anywhere fix. Pattern matching is necessary starting at −145 dbm and can be carried on down to as low as −150 dbm. A z-count or intMs can be thereby obtained so the SV can be used in an anywhere fix. However at such signal levels, the ephemeris needs to be obtained over a network from a server, or alternative source. Below signal levels of −150 dbm, the NAV-data is not reliable enough for a pattern match. The NAV-data must be obtained from the server 204, and the SV's with such weak signals can only participate in a fix when the uncertainty is less than 150-km.

During initial SV acquisition, ephemeris level accuracy is not needed. An almanac, or downgraded ephemeris, is adequate to predict the data needed for preposition. Ephemeris level accuracy is also not needed for fixing. A time-out is defined for the ephemeris age for positioning. Such threshold can be relaxed and still maintain respectable positioning if the accuracy degradation as a function of time is properly modeled. The age threshold can be a controllable parameter so that the customer can select the desired level of performance.

NAV-data subframe data from the server is needed for first fix or setting time. After that, subFrames are no longer requested by a client. The NAV-data decoded by the client can be sent to the server for the server to do the pattern-matching.

A server connection is not needed by an OMNI client when there are three or more SV's with signal levels all better than −145 dbm. The time-to-first-fix (TTFF) will be longer if the ephemeris must be collected. In some cases, previously collected ephemeredes can be used.

A server connection is not needed by an OMNI client when previously collected ephemeredes for SV's are on-hand and sigmaPos is less than 150-km. The minimum numbers of SV's needed depends on sigmaTime. Such time uncertainty can be reduced with a real-time clock (RTC) that is software-compensated for temperature drift. So three SV's are needed with such RTC, and four SV's without the RTC

Otherwise, solving for a fix will require that the OMNI client contact the server and request certain information. The NAV-data subframes are needed when the SV's signals are −145 dbm to −150 dbm and sigmaPos>150-km. The intMs for these SV's is needed for them to participate in the first fix. If only three −145 dbm or weaker SV's are available and no other better means of accurate time, pattern matching can be used. A so-called no-Z fix with four SV's is then used.

The ephemeris has to be requested when SV's signals are not stronger than −145 dbm and their ephemeris has timed out. In such case, the fastest TTFF possible is desired.

A main program application can periodically turn-on the GPS receiver and get a fix. Such decides how far the receiver has moved since the last fix, or decides simply if the GPS receiver has left a predefined zone. The time interval between fixes is selected to keep sigmaPos within 150-km so that the intMs is not needed on weak SV's under −145 dbm. Such extends the ability to retain high sensitivity fixes without needing a server connection to request NAV-data subFrames. The timing of the server requests is adaptive. This is needed to provide a quiet client/server connection when there is adequate performance without it.

The OMNI client must evaluate the data it has, the data's age, and the likelihood of acquisition success, e.g., number of SV's and the signal levels. The OMNI client then decides whether to make the connection and what data to request. The adaptivity can be disabled and the server connection can made by explicit commands. A master application may decide to make a server connection every one hour. Thus, for fixes done every five minutes, the twelfth one will make a server connection.

A broadcast type ephemeris service can be used where the master application collects the data and then pushes it into the client through generalized API's. The client can be empowered to make a server connection any time during a session.

The client typically has a time uncertainty larger than one millisecond before a fix. So a conventional search of the complete codephase is done prior to the first fix. After the first fix, or with some restarts using RTC, the client may have a time uncertainty under one millisecond. In these cases the client can reduce the code search window. So even if the client has time from a measured z-count, the client still has to search all 1023-chips of code uncertainty.

Critically, in all embodiments of the present invention, a portion of DSP 216 (FIG. 2) is kept alive after power-down and is maintained as the time source between client sessions. During power-down, power is actually turned off to RF ASIC 212 ASIC 212, sampleClock, sampleMem, and OSM. But the power is kept alive to the crystal oscillator 218 and a millisecond-interrupt internal to DSP 216 so it can process it.

The client 202 has the advantage of the oscillator stability in maintaining the time calculated in a last position fix for a relatively long period until the next GPS session.

The client empirically determines the accuracy that can be achieved. To do this, the client can set time from a fix, and then stop tracking for a given off time. Then the client re-enable tracking and use the measured z-count as an indication of how far the clock has moved. The clock adjustments the client perform in the fix will tell us how far the clock has moved with respect to GPS time.

This technique requires that the client maintain power to the crystal oscillator 218 and DSP 216. It requires that the client has the ability to keep power to the GPS circuits even when the device is apparently turned off.

In a typical shutdown, when all firmware state machines are idled, all “int1” activity should cease. The client will also disable automatic output so the “int0” activity will also stop. At this point, there should not be even any background activity and the DSP 216 should finish its tasks and fall to halt mode. The halt mode is a special low power feature of the DSP 216 where it remains idle until there is an interrupt at which time it automatically jumps to that interrupt service routine.

When the millisecond-interrupt arrives, the DSP 216 wakes up and services it, advancing the all timers, including GPS time estimates, by one millisecond. The DSP 216 will then go back to halt mode when it exits the interrupt.

The client needs to output absolute time. So whenever the GPS-TOW (in seconds) rolls from 604799 to 0, the client should also increment a week counter. The client currently don't maintain absolute GPS week on DSP 216 so that feature must be added for the RTC mode. None of the processing on DSP 216 uses the week number however. It would only be included to provide absolute GPS time outside of DSP 216. The client will also benefit by maintaining a time-off variable that can be read by the client when it restarts.

The cycle is completed when the another client session commences. At this point, the client will request continuous time message output again and the client will be able to observe whether DSP 216 has remained on or off, and whether it has confidence in the time. Thus, the client will need to construct a bias and drift uncertainty that can be used by the client to decide how it can use the time from DSP 216.

When the client sends a time adjustment, it should also send the estimated accuracy of the time set. This accuracy is accepted by DSP 216. This accuracy will be automatically propagated inside DSP 216 according to a model of the stability of the crystal with respect to temperature. If there is more temperature change, the bias and drift uncertainty should propagated at a higher rate. However, if the temperature is constant, the uncertainty is propagated at a slower rate.

When the client is finishing a session, it will communicate to DSP 216 to disable all the GPS tracking functions. This will shut off RF ASIC 212, and disable all GPS state machines. The client can select either deep-sleep mode or RTC mode. If deep-sleep is selected, then the client turns off the crystal oscillator 218. However, if RTC mode is requested, then the client keeps on the crystal.

The client can also request that DSP 216 maintains its millisecond clock. It can also decide to shut off everything, including the millisecond if it desires to produce the lowest power consumption. Most power is consumed only when transistors change state, and thus, if all circuits are off, the consumption is minimized. In deep sleep, the oscillator is also turned off when the DSP 216 detects deepSleep mode and encounters the halt command.

If the client turns off DSP 216 but requests the millisecond option, then the firmware disables all modes that it can, but keeps the millisecond and oscillator active. The client also disables the automatic output of time messages to keep DSP 216 as quiet as possible.

In embedded applications there is always the possibility that the power will be turned off without using a proper shutdown procedure. Ideally, the client would get a turn-off request from the application using the client. This starts the chain of events in the following table.

The client receives a request from host application that it wishes to turn off GPS; The client sends messages to DSP 216 to turn off all GPS tracking engines            (ODSM, IDSM, TSM); The client sends the preferred shut down mode of DSP 216 and waits for            acknowledgement of the action. If it wants the lowest power mode, it           requests DEEP_SLEEP_MODE. If it wants the most information, it           requests RTC_MODE; The client flushes all data from its serial buffers. The client finishes processing all events, and thus, the key memory elements that are            battery backed become static; The battery backed data is saved; The client sends a message back to host application that it is in a safe state to shut down;      and The host de-allocates the client applications.

Such procedure is a preferred method when the shutdown occurs in a controlled way.

It is also possible that a less orderly shutdown will occur. In such case, the host application abruptly terminates the client code before it can perform its shut down procedure. E.g., the DSP 216 was not properly shutdown, or the power was removed from the system. This could occur when someone removes the battery or the battery cannot provide enough current to maintain proper operating levels.

If the DSP 216 is still being powered, a checking logic can be included to detect if communication with the client has terminated. If so, the DSP 216 does its own power-down sequence. If the DSP 216 was simply powered down, when it turns on again it will be in its default state and not a problem. However, there is a risk that the client was writing to a section of the battery-backed file at the time and the data may be corrupted. A variable could be written to flag to indicate the DSP 216 shut down properly. If the variable is not set when the data is read at startup, then the data is considered invalid and is set explicitly to an invalid state.

During turn-on, the client application is started and tickles DSP 216 to wake it up and requests the time-status. Such tickle also immediately starts the TCO in DSP 216. When the time-status is received, the client sets its bias and drift uncertainty to unknown states if the time-status is invalid. In this way, any battery-backed data or data from the server can properly battle themselves for the smallest sigma. The GPS time will also be zero at this point until the client interrogates its other information sources.

If the time-status is valid, then the client copies the time uncertainty data and initializes its internal data structures with these data. If better data is available from other sources, then it can battle those uncertainties when that data becomes available. The client then reads the status concerning the battery backed data. It assimilates data sections that have the data validity flags set. Then the client is ready to evaluate its available information and begin requesting data from the server.

During an orderly shutdown, if RTC_MODE is selected, the time data is sent to DSP 216. It waits for acknowledgement as part of the shutdown procedure. If the time data is only sent in the shutdown procedure, a non-orderly shutdown of the client occurs where the power remains ON to DSP 216. In this case, the RTC_MODE will not be enabled, and DSP 216 continues to run. The state machines will still be running, and the time information will be lost.

A periodic handshake is preferred between the client and DSP 216. If that handshake is lost, and DSP 216 is still powered, then DSP 216 needs to put itself in a quiet state.

The client could still keep time active if the timing data can be sent periodically. However, maintenance of the timing data will be expensive computationally so the client will need to orchestrate these computations to be a very low priority.

The following time and uncertainty maintenance equations can be used in designs. At some point, the client will get time information from either a position fix, a z-count, a pattern-match, or from another source such as a different real-time-clock, or from the server latency estimator. The client can generate a model of the accuracy of each data source for both bias and drift. The calculations are done meters for bias, and m/s for drift. An SCXO model used has units of m/s. In prepositioning, the client always computes in pseudorange, and pseudorange rate, and then convert to hardware units of 64^(th) chips and carrier NCO units. So, if the client wants to use the number directly in DSP 216, then the client will have to convert to NCO units.

$\begin{matrix} {{{{{Drift}\left( {m\text{/}s} \right)}/{\lambda \left( {m\text{/}{cycle}} \right)}}*{{bits}/{Hz}}} = \begin{matrix} {{drift}*\left( {1575.42\; e\; {6/2.99792458}\; e\; 8} \right)*} \\ {{2^{24}/528}\text{,}000^{\;}} \end{matrix}} \\ {= {{drift}*166.9789113}} \end{matrix}$ Drift  NCO  units  (bits/Hz = drift(m/s) * 166.9789113

The client should periodically send the following atomic data set:

1. The current bias (b_(fix)) in units of meters with an LSB of 2⁻⁶ 2. The current drift (d_(fix)) in units of meters/sec with an LSB of 2⁻¹² 3. The bias sigma (σ_(bf)) in units of meters with an LSB of 2-6 4. The millisecond sigma (σ_(millisecond)) in units of milliseconds with an LSB of one millisecond 5. The drift sigma (σ_(df)) in units of meters/sec with an LSB of 2-12 6. The (TCO_(fix)) measurement closest to estimate of the drift (TCO_(fix)) in units of counts with an LSB of 7. The msec16 that corresponds to the bias exactly and also all the other parameters

For the frequency uncertainty, the client has just one parameter σ_(df). However, for the time uncertainty, the client prefers to work with two parameters. In many cases, before the first fix, the time uncertainty is many milliseconds, so the client prefers to break the time uncertainty into two components: the part that is multiple millisecond, and the part that is sub-millisecond.

Accuracy of bias and drift from different sources

Fix type or information The millisecond sigma The bias sigma (σ_(bf)) The drift sigma (σ_(df)) source (σ_(millisecond)) (milliseconds) (meters) (m/s) From DSP 216 RTC Computed from data Computed from data Computed from data sent to DSP 216 sent to DSP 216 sent to DSP 216 originally originally originally Approximate time from 1000 300,000 Based on SCXO RTC or client device availability or from battery backed drift estimate Approximate time from 1000 300,000 Based on SCXO server availability or from battery backed drift estimate Accurate time from SigmaServer (calculated 300,000 Based on SCXO server in the server latency availability or from logic) battery backed drift Probably best case 10-20 msec. estimate Time from a pattern 10+. (There is a small 300,000 Based on SCXO match probability this time availability or from could be wrong) battery backed drift estimate Time from a z-count  10 300,000 Based on SCXO availability or from battery backed drift estimate Time from a no-z fix >> 10. Could be as large TDOP * σ_(measPR) TDOP * σ_(measPRR) as 500, but is a function Paradox state: bias With good of geometry and meas sigma small, measurement, << noise millisecond sigma 10 m/s (indoor) larger << 1 m/s (outdoor) Time from fix that   0 TDOP * σ_(measPR) TDOP * σ_(measPRR) solves for time error Paradox state: bias With good linearized around time sigma small, measurement, << from pattern match, or millisecond sigma 10 m/s (indoor) z-count, or DSP 216 larger << 1 m/s (outdoor) RTC mode or confident server time

The propagation equations update the drift with the TCO and previous drift, propagate the bias with the drift, do millisecond adjustments when the bias exceeds +/−½ millisecond, update the drift uncertainty, and update the bias uncertainty. In order to propagate the bias and bias sigma, the client will need to update the drift and drift sigma. Fundamentally, the bias is the integral of the drift, and the bias sigma is the integral of the drift sigma.

The client estimates the drift based on the drift data at the fix and the current SCXO data based on current TCO reading.

driftHat = function(drift  at  fix, SCXO  drift  at  fix, current  SCXO  drift)  

The bias with a drift estimate is propagated. When the bias exceeds +/−½ msec, the client adjusts the DSP 216 millisecond so that the millisecond stays accurate

Bias = bias + driftHat * dt If (bias > millisecond/2) {    Bias −= millisecond    Adjust DSP 216 millisecond +1 millisecond } else if (bias < −millisecond/2) {    bias += millisecond    Adjust DSP 216 millisecond −one millisecond }

The client updates the drift sigma based on the drift at the fix and the current drift model based on current TCO.

driftSigma(t) = function(drift  sigma  at  fix, SCXO  sigma  at  fix, current  SCXO  sigma) 

The client propagates the bias uncertainty as the integral of the drift uncertainty.

biasSigma = biasSigmaAtFix + integral(driftSigma(t)) * t

When the next session begins, the client sends back the bias, drift, and sigmas. Hopefully, the bias sigma will grow to less 10 or 20 millisecond meaning the client can actually do a fix without reliance on time from any other source.

Thus, the following parameters are defined from the observations:

1. δ(t_(fix))=d_(fix)−d_(scxo)(t_(fix)). This signed parameter is measured at fix time and represents the model error at the fix time. It remains constant until the process is re-started.

2. δ(t_(now))=d_(scxo)(t_(now))−d_(scxo)(t_(fix)). This is a signed parameter and represents the change in the model from fix time.

The first goal is to propagate the bias with the best estimate of the drift. Of course the client also has the requirement that the estimated bias sigma is the least upper bound (LUB) of the true bias error. In this way, as the client propagates the bias sigma it will remain sufficiently small so that after a long period of time the sigma indicates that the client can use the bias to compute a fix. Stated conversely, if the bias sigma is too large, then the client cannot confidently use the bias estimate to compute a fix.

Around the time of the fix the measured drift is certainly the best estimate the client has. As long as the temperature stays near the fix temperature then the client can use the fix drift as the propagation source. This will also produce the smallest drift sigma, and thus, the bias sigma will integrate slowly and help produce the LUB criteria.

Conversely, when the time since the last fix is large, the client tends to lose confidence in the fix drift as a reference point and the client prefers to transition to the completely modeled data. However, the client has to use the larger sigma to propagate.

In between these two extremes, it seems logical to blend the two sources (fix and model) together. Empirically, the client has seen that it is common for the drift trajectory to be somewhat biased from the SCXO model for long periods. However, from day to day, it is equally likely that the bias is above or below the model. If the client blends the data together, and leans on the fix drift as long as possible to produce the current days bias, then the client can hopefully obtain a more accurate drift and also reflect that in a smaller drift sigma.

One key observable is |δ(t_(fix))|, the difference between the measured and modeled drift at the fix time. Assume that the client will perform some integrity checking on the measurements used in the fix, and also the fix itself. For example, if the computed velocity is close to zero, this is like a stable point and it is likely that the Doppler measurement errors are small. Thus, the client will assume that the measured drift from the fix has high confidence. Degradations due to using indoor measurements will be reflecting the drift sigma from the fix. From this point of view, the client tends to have more confidence in the recent measured data.

The aging or shock of the crystal will tend to produce a bias between the crystal and the model. The client must accept there will be discrepancies between the measured and modeled data. The client will develop a methodology to propagate the bias and bias sigma with a drift and drift sigma that has a stronger weighting to the measured data but slowly converges towards the modeled drift as the time from the measurement increases.

A simple formulation is defined for migrating between the drift estimated at the fix to a drift estimate from the SCXO. The client will control the rate of the migration via a couple simple parameters.

The client will be blending two estimates of drift, one that is the SCXO model based on the current TCO, and the other is the adjusted fix drift. The adjusted fix drift is simply the drift from the fix corrected for the change in the model due from the current time and the time of the fix.

The adjusted fix drift is:

d _(a)(t _(k))=d(t _(fix))+d _(scxo)(t _(k))−d _(scxo)(t _(fix))

The client defines the blended drift estimate as:

d̂(t _(k))=[1−α]*d _(a)(t _(k))+α*d _(scxo)(t _(k))

Notice that the client has two estimates of drift that the client blends with the parameter α. This parameter provides,

1. When t_(k)=t_(fix), then the client want d̂(t_(k))=d_(a)(t_(k)). This requires α=0.

2. When t_(k)>>t_(fix), then the client want d̂(t_(k))=d_(scxo)(t_(k)). This requires α=1

Thus, defining the parameter N, the count of the times the filter has been updated,

Define:

N₁=1, if k=1

N_(k)=k if k<=N_(max)

N_(k)=N_(max) if k>N_(max)

The client also wants to be able to perform the multiplies using integer math, thus, the client defines,

α=N _(k)/2^(M)

1−α=(2^(M) −N _(k))/2^(M)

N_(max)=2^(M)

The client only needs to specify one parameter M and then the rate of change is completely defined.

M N_(max) = Seconds to final value Hours to final value 8 256 .071 10 1025 .284 12 4096 1.14 14 16368 4.55 16 65536 18.2

The same formulation is used for the drift sigma.

σ_(d)(t _(k))=[1−α]*σ_(fix)+α*σ_(scxo)(t _(k))

The models are updated every second, so there is an implicit multiply by one second, the client can update the bias as follows:

b(t _(k))=b(t _(k-1))+d̂(t _(k))

The client will also update the bias sigma as the integral of the drift sigma, so there is an implicit multiply by one second as follows:

σ_(b)(t ₁)=σ_(bfix)

σ_(b)(t _(k))=σ_(b)(t _(k-1))+σ_(d)(t _(k))

The crystal oscillator 218 and DSP 216 hardware can be used to verify that accuracy of the model and also select the correct time constant. Real DSP 216 hardware tracks at least one satellite so that the client can always measure true bias and drift. The client does this by linearizing the pseudorange and pseudorange rate around the true position as follows:

a. LPR=intms*cmsec+codePhase−rangeHat−corr

b. LPRR=rangeRate−rangeRateHat

c. The model is that LPR=biasTrue and LPRR=driftTrue

The very first measurement, the client sets the receiver clock terms to null out LPR and LPRR. Any component of LPR larger than +/−½ millisecond goes into millisecond clock adjustments. The sub-millisecond portion residual is the starting bias. After the first adjustment, the client now starts to form bias and drift from the SCXO data and equations given above. If the bias grows over the +/−½ msec range, the client does the clock adjustment as usual. The LPR and LPRR will provide the true bias and drift and the equations will provide the estimated bias and drift. Also, when ever there is measured z-count, if the client corrects this with the true intms rather than assuming 70, the client will have the true millisecond time error.

FIG. 3 illustrates a multiple re-use of frequencies in an integrated combination of a satellite navigation receiver and communication device system embodiment of the present invention, referred to by the general reference numeral 300. The system 300 has a GPS part 302 and a cellphone part 304 that are both integrated within the same portable or handheld device. The GPS part 302 includes a GPS antenna 306, a GPS RF-receiver 308, a GPS digital processor and frequency generator (chip) 310, a GPS oscillator (oscGPS) 312, and a low frequency oscillator (osc3) 314. A clock selector 316 selects between normal and low-power operation, e.g., by clock frequency choices. The cellphone part 304 includes a cellphone antenna 318, a communications transceiver 320, a digital signal processor 322, a phone CPU-power control-clock selector 324, a divider 325, and a voltage controlled oscillator (VCO) 326.

The divider 325 is included to increase the relative stability and accuracy of the reference frequency applied to VCO 326. The cellphone part 304 is not as dependent on its cellphone base station providing good reference frequencies for proper CDMA operation when it has very accurate references locally. If all the mobile phones operating in a base station area had such GPS-provided reference signals, the base station would not need to be equipped with expensive atomic clocks.

Communication devices generally need their own faster processor crystals, since the typical VCO oscillator operates in the relatively slow 10-13 MHz range. Some processors need to be able to switch to even higher frequencies for special communications functions. An intermediate frequency (IF) used in the down conversion process and produced in the RF chip can be used to run applications that require more computation resources. A low frequency oscillator is selected for low-power time maintenance, and is typically a 32-KHz crystal.

At turn-on, the GPS chip 310 loads a default startup processor clock selection for the phone CPU 324, e.g., from a non-volatile memory location. The host processor clock frequency is generated by multiplying GPS clock for input to a numeric controlled oscillator (NCO). Such NCO can then be digitally programmed to generate any output frequency up to half of its input frequency. The GPS part 302 puts itself into a sleeping, low-power mode and lets the user decide when to use the GPS. A host system begins operating when its clock is stable. The host can then chose different frequencies, via the communication between the host CPU 324 and GPS digital chip 310. The phone part 304 can request a general frequency number by sending a desired frequency. It can control when the VCO is turned onto the frequency to be used, e.g., 12.6 MHz for the Japanese PDC system. If the GPS part 302 is tracking GPS satellites and solving for its frequency error, it can then compensate the requested 12.6 MHz by what it knows to be the error in the GPS crystal. Thus, it can provide a more stable frequency to the communications device VCO 326. If low-power mode is needed, the communications device 304 can request that the GPS part 302 enter low-power mode. It keeps time with the low frequency crystal input from osc3 314. If the GPS oscillator 312 can be made low-power, the 32 kHz oscillator 314 can be eliminated. The GPS part 302 can send interrupts to the phone part 304 on time-events lines to wake certain processes that need to occur at regular but accurate intervals.

The GPS part 302 provides a VCO frequency corrected by GPS satellite time standards, and is very accurate. A variable VCO from the GPS part 302 can be supplied that is responsive to requests from the phone part 304, e.g., for frequency stepping operations needed in frequency-division multiple access (FDMA) systems. An analog sinusoidal voltage can be output by the GPS part 302 to approximate a crystal output. Such may be requested by using bits from the top of an adder for accurate phase. Such phase is converted with a table having more representative levels for a sinusoid rather than the linear counter value. A table value can be converted to an analog signal in a digital to analog converter. The overflow of the adder can also be used to generate a simple 1-bit digital clock frequency. Time interval interrupts can be constructed from combinations of the GPS second and millisecond interrupts in the GPS receiver, and osc3 314. The timing of events can be phased to any time frame of reference by using offsets.

In a method embodiment of the present invention, the number of source oscillators in an integrated combination navigation receiver and cellphone is reduced to two, e.g., a GPS crystal oscillator at about 27-MHz and a watch crystal oscillator at about 32-KHz. A multiplier is connected to the GPS crystal oscillator to produce higher frequencies. Two numeric controlled oscillators (NCO1 and NCO2) are used, respectively, to generate VCO and host CPU frequencies. A time event logic produces time events to the host CPU from combinations of GPS msec interrupts, GPS second pulses, digital offsets, and the watch crystal oscillator.

In general, embodiments of the present invention improve both manufacturing costs and device performance. For example, the navigation receiver supplies a reference crystal frequency to the communications device to re-use a crystal, thus eliminating the need for a second crystal, the communications VCO crystal. And, when the navigation receiver supplies such reference crystal frequency to the communications device, the communications receiver sensitivity is improved because the frequency uncertainty is so much reduced the initial frequency search space can be trimmed. The receiver is thus able to search for signal in the frequency domain using smaller increments or steps, but still be able to find initial lock in a reasonable time.

Although the present invention has been described in terms of the presently preferred embodiments, it is to be understood that the disclosure is not to be interpreted as limiting. Various alterations and modifications will no doubt become apparent to those skilled in the art after having read the above disclosure. Accordingly, it is intended that the appended claims be interpreted as covering all alterations and modifications as fall within the “true” spirit and scope of the invention. 

1. A method for limiting the number of crystals needed for reference oscillators when integrally combining a navigation receiver and cellphone to two, comprising: providing a GPS oscillator with a first crystal operating at about 27-MHz; providing a watch oscillator with a second crystal operating at about 32-KHz; connecting a multiplier to the GPS oscillator to produce higher frequencies to operate a GPS receiver; generating a VCO frequency with a first numeric controlled oscillator (NCO1) to operate a communications transceiver; generating a host CPU frequency with a second numeric controlled oscillator (NCO2); and producing time events to host CPU with logic combinations of GPS-msec interrupts, GPS-second pulses, digital offsets, and said watch oscillator. 